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Abstract 

Estimation of the minimum eigenvalue of a quantum Hamiltonian can be formalised as 
the Local Hamiltonian problem. We study the natural special case of the Local Hamiltonian 
problem where the same 2-local interaction, with differing weights, is applied across each pair 
of qubits. First we consider antiferromagnetic/ferromagnetic interactions, where the weights 
of the terms in the Hamiltonian are restricted to all be of the same sign. We show that for 
symmetric 2-local interactions with no 1-local part, the problem is either QMA-complete 
or in StoqMA. In particular the antiferromagnetic Heisenberg and antiferromagnetic XY 
interactions are shown to be QMA-complete. We also prove StoqMA-completeness of the 
antiferromagnetic transverse field Ising model. Second, we study the Local Hamiltonian 
problem under the restriction that the interaction terms can only be chosen to lie on a 
particular graph. We prove that nearly all of the QMA-complete 2-local interactions remain 
QMA-complete when restricted to a 2D square lattice. Finally we consider both restrictions 
at the same time and discover that, with the exception of the antiferromagnetic Heisenberg 
interaction, all of the interactions which are QMA-complete with positive coefficients remain 
QMA-complete when restricted to a 2D triangular lattice. 


1 Introduction 

Calculation of the ground-state energy of quantum Hamiltonians with interactions obeying 
locality constraints is a fundamental problem in physics. This is encapsulated within quantum 
information theory as the LOCAL HAMILTONIAN problem. A Hamiltonian H on n qubits (a 
Hermitian matrix acting on is said to be fc-local if it can be written as H = 

where each interaction term Hj acts non-trivially on at most k qubits. Let be the ground- 
state energy (the lowest eigenvalue) of H, and assume that ||iL; || ^ poly(n) for all j. Then the 
Local Hamiltonian problem is to determine whether X{H) ^ a, or X{H) ^ b, for some a, b 
such that b — a ^ 1/ poly(n). 

It was shown by Kitaev [17] that, even for k = 5, Local Hamiltonian is QMA-complete. 
QMA is a complexity class which is the quantum analogue of the classical complexity class 
NP [24]. A problem is said to be “complete” for a complexity class if it is contained within 
that class, and any other problem in that class can efficiently be reduced to it; if a problem 
is QMA-complete, it is likely that there is no efficient quantum (or classical) algorithm to 
solve it. Kempe, Kitaev and Regev [16] later showed that the Local Hamiltonian problem 
remains QMA-complete even for k = 2, whereas the special case where k = 1 can easily be 
solved efficiently. The Hamiltonians occurring in this hardness proof are somewhat artificial. 
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A number of subsequent works have therefore attempted to show that Local Hamiltonian 
remains QMA-complete, even with more physically realistic restrictions on the interactions [22, 
3, 28, 27, 12, 10]. This simultaneously makes the theory more useful for applications and 
provides insight into the crucial features of physical theories that determine their complexity. 
For example, it allows one to ask, and hopefully answer, questions of the form “is the Heisenberg 
model on a square lattice more complex than the Ising model on a triangular lattice?”. 

A particularly natural restriction to consider is to restrict the types of interactions allowed. 
Let 5 be a set of allowed interaction terms on at most 2 qubits. The 5-Hamiltonian problem 
is a restriction of LOCAL HAMILTONIAN where k = 2 and the Hamiltonian H = ^2 OiiHi can be 
written as a linear combination of terms Hi € S, and where each ai is a positive or negative real 
weight such that \ai\ ^ poly(n). For example, the general Ising model, where H = Ylij 
corresponds to S = {ZZ}. Similarly, the (general) Heisenberg model is {XX + YY + ZZ}] the 
XY model is {XX + YY}; the Ising model with transverse magnetic field is {ZZ,X}. Note 
that X, Y, Z are the Pauli matrices and we omit tensor product symbols for readability. 

A complexity classification of the 5-Hamiltonian problem was obtained in [12]. This can 
be stated as follows: 

Theorem 1 (Cubitt and Montanaro [12], Bravyi and Hastings [8]). Let S be an arbitrary fixed 
subset of Hermitian matrices on at most 2 qubits. Then: 

• If every matrix in S is 1-local, 5-Hamiltonian is in P; 

• Otherwise, if there exists U G SU{2) such that U diagonalises all 1-qubit matrices in S, 
and diagonalises all 2-qubit matrices in S, then 5-Hamiltonian is NP-complete; 

• Otherwise, if there exists U G SU{2) such that, for each 2-qubit matrix Hi G S, 

aiZ®"^ -\- Ail -\-1 Hi, where a* G M and Ai, Hi are arbitrary single-qubit Hermitian matrices, 
then 5-Hamiltonian is StoqMA-complete; 

• Otherwise, 5-Hamiltonian is QMA-complete. 

The third case in Theorem 1 was originally shown to be complete for a complexity class corre¬ 
sponding to the transverse Ising model in [12]; this was later sharpened to StoqMA-completeness 
by Bravyi and Hastings [8]. The Local Hamiltonian problem can be seen as a quantum 
generalisation of classical constraint satisfaction problems, with each term in the Hamiltonian 
corresponding to a constraint. From this perspective, Theorem 1 is a quantum analogue of a 
classical dichotomy theorem of Schaefer [26] classifying the complexity of constraint satisfaction 
problems in terms of the types of allowed constraints. 

Each of the above complexity classes (P, NP, StoqMA, QMA) corresponds to a computa¬ 
tional model. P is polynomial-time classical computation (decision problems that can be solved 
efficiently by a classical computer). NP (“nondeterministic polynomial-time”) corresponds to 
decision problems whose solutions can be checked efficiently. StoqMA is a complexity class intro¬ 
duced by Bravyi, Bessen and Terhal [6] whose definition as a computational model is somewhat 
technical, but which corresponds to the special case of the Local Hamiltonian problem where 
the matrices Hj have non-positive off-diagonal elements (known as “stoquastic”). Finally, QMA 
is the complexity class corresponding to decision problems whose solutions can be checked in 
polynomial time by a quantum computer. We have P C NP C StoqMA C QMA and it is 
believed that each of these inclusions is strict. Some natural examples of interactions falling 
into each of these classes: the Ising model is NP-complete [1]; the transverse Ising model is 
StoqMA-complete [8]; and the Heisenberg and XY models are QMA-complete [12]. 
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Although Theorem 1 gives a precise and complete classification of the complexity of 2-qubit 
interactions, it suffers from several shortcomings which make it less meaningful as a statement 
about physics: 

1. The interactions do not have any constraints on their spatial locality. They can occur 
across large distances and each individual qubit can interact with arbitrarily many others; 

2. The interaction terms Hj can appear with either positive or negative signs in the final 
Hamiltonian H, corresponding to two quite different physical interpretations; 

3. The weights of individual interactions can be as large as poly(n), and the precision with 
which we are asked to estimate the ground-state energy scales as 0(1/ poly(n)). In physical 
systems the interactions are usually of weight 0 ( 1 ). 

In this work we address the first and second of these issues. We do not address the third (but 
see [9, 10, 11] for some cases where QMA-completeness has been proven with interactions of 
strength 0 ( 1 )). 

1.1 Statement of results 

We study the 5-Hamiltonian problem with physically motivated restrictions. The first re¬ 
striction is the additional constraint that all weights a* must satisfy ai ^ 0. We will label this 
problem 5’^-Hamiltonian, although for the majority of this paper we consider sets of only one 
element. Physically, this restriction allows us to characterise the complexity of models which 
have only either antiferromagnetic or ferromagnetic interactions. This restriction is also useful 
when viewing the 5-Hamiltonian problem as the quantum analogue of a weighted constraint 
satisfaction problem. Indeed, from this point of view it is difficult to interpret the meaning of 
negative coefficients. 

Our main result in this setting is the following theorem: 

Theorem 2. Given a 2-qubit interaction H = aXX -t- /3YY 'yZZ, the problem 
Hamiltonian is either: 

i) QMA-complete, if a (5 > 0, a-|- 7>0 and /3 -t- 7 > 0; 

ii) in StoqMA, otherwise. 

Furthermore, if a = —f3 ^ 0, a -|- 7 > 0 and (3 -\- j > Q , then {//}'''-Hamiltonian is StoqMA- 
complete. If a = j3 and 7 ^ —|q:|, {f^}■'■-HAMILTONIAN is in P. 

The latter part of the theorem can easily be seen to hold under relabelling a, (3, 7 . It 
was shown in [ 12 ] that if at least two of a, (3, 7 are non-zero, then the corresponding problem 
{H}-hamiltonian is QMA-complete; and so Theorem 2 shows that the restriction of having 
positive coefficients does indeed make a significant difference to the complexity. 

Observe that we can transform any 2-qubit interaction H which is symmetric under swapping 
the qubits, and has no 1-local part, into the form required for Theorem 2 using local unitaries 
(see [12] and references therein). Theorem 2 thus classifies the complexity of all symmetric 
2-qubit interactions with no 1-local part. We illustrate this classification in Figure 1. 

Important special cases of Theorem 2 are the antiferromagnetic Heisenberg model {H = 
XX -\-YY ZZ) and the antiferromagnetic XY model (H = XX -\-YY), which are both shown 
here to be QMA-complete for the first time. A similar result has been shown very recently 
for the antiferromagnetic XY interaction by Childs, Cosset and Webb [11], using different 
techniques. Their result holds even when the a* coefficients are all either 0 or 1, but only at 
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Figure 1; The complexity of interactions of the form H = aXX + I3YY + ^ZZ. The black 
line indicates StoqMA-complete interactions. The origin is in P in the first two cases, and 
NP-complete in the third. 


fixed magnetisation (i.e. within a fixed eigenspace of Z®”). Another interesting case which 
does not fit into the framework of Theorem 2 is the antiferromagnetic transverse Ising model, 
S = {ZZ, A}. We show that this case of 5'’“-Hamiltonian is also StoqMA-complete. 

The second restriction we study is to require that all interactions only occur on the edges 
of a particular interaction graph, where each vertex corresponds to a qubit. We focus on two 
natural graphs: 2D square and triangular lattices. Square lattices have been studied by several 
authors. It was shown by Oliveira and Terhal that general 2-local qubit Hamiltonians on a 
square lattice are QMA-complete [22]. Then the Heisenberg interaction with arbitrary 1-local 
interactions was shown to be QMA-complete on a square lattice by Schuch and Verstraete [28]. 
It was shown in [12] that any QMA-complete 2-qubit interaction remains QMA-complete on the 
square lattice, if arbitrary 1-local terms are also allowed. However, no interaction has previously 
been shown to be QMA-complete on a square lattice without 1-local terms. 

Our first result in this setting concerns square lattices without sign restrictions. To state 
the result, it will be helpful to decompose a general 2-qubit interaction H into Pauli matrices: 
H = Mijai (8> crj+ 1-local terms. Then the Pauli rank of H is the rank of the 3x3 matrix 

M. 

Theorem 3. Let H be a 2-local qubit interaction with Pauli rank at least 2, such that the 2-local 
part of H is not proportional to XX -\-YY -\-ZZ. Then {Pf}-HAMiLTONiAN is QMA-complete, 
even when the interactions are restricted to a 2D square lattice. 

We believe that the one excluded case in this theorem, the Heisenberg interaction, is in 
fact QMA-complete; however, we have not been able to prove this. QMA-completeness of 
the Heisenberg interaction would follow if there existed an exactly solvable special case of the 
Heisenberg model satisfying certain constraints (see below). If this final case were found to 
be QMA-complete, we would have shown that all of the QMA-complete 2-qubit interactions 
remain QMA-complete on a square lattice. 

Our second result about interaction geometry concerns 2D triangular lattices with sign 
restrictions. A square lattice can be viewed as a sublattice of a triangular lattice, and so any 
hardness results on a square lattice such as Theorem 3 can be trivially extended to the triangular 
lattice. However, we can obtain stronger results by working with triangular lattices directly. 

Theorem 4. Let H = aXX j3YY yZZ he a 2-qubit interaction such that a -|- /3 > 0, 
q:-|-7>0,/3-|-7>0 and H is not proportional to XX-\-YY -\-ZZ. Then {Hj^-HAMILTONIAN 
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is QMA-complete, even if the interactions are restricted to the edges of a 2D triangular lattice. 

An important special case of Theorem 4 is the antiferromagnetic XY model on a triangular 
lattice, which we show is QMA-complete. Observe that Theorems 2 and 4 overlap, but Theorem 
4 is not quite a stronger result than Theorem 2. We do not expect to be able to show QMA- 
completeness results for {i?}"*’-HAMiLTONiAN for any of the interactions H in Theorem 4 when 
restricted to a square lattice. Indeed, for any bipartite interaction graph, we can perform a 
local change of basis conjugating by Z on every qubit on one side of the bipartition. This will 
effectively send a, P ^ —a, —/?, thereby converting any of the QMA-complete Hamiltonians 
from Theorem 2 into one for which approximating the ground state energy is contained in 
StoqMA. Together, Theorems 3 and 4 thus highlight a way in which the triangular lattice can 
rigorously be said to be more complex than the square lattice. 

1.2 Proof techniques 

The main proof technique used is to apply perturbative gadgets [16, 22, 8, 12] to show re¬ 
ductions from problems that are already known to be QMA-complete. The basic idea is to 
use a Hamiltonian H = AHq + V, for some large A, such that the low energy sector of H is 
approximately equal to the ground space of Hq, and that within this space the action of H is 
determined by V to some order in a perturbative expansion in 1/A. This allows new interac¬ 
tions to be approximately implemented that were not previously available. For Theorem 2, a 
second order perturbation “sign” gadget is used to allow a Hamiltonian with positively weighted 
terms to simulate a Hamiltonian with mixed signs that is known to be QMA-complete. The 
StoqMA-completeness results are based on perturbative reductions from the transverse Ising 
model, which was recently shown to be complete for StoqMA by Bravyi and Hastings [8]. 

The main technical challenge in order to prove QMA-hardness for interactions on a 2D 
lattice turns out to be proving QMA-hardness of the XY model. This proof consists of two 
steps. It is known that Hamiltonians made up of XX and ZZ interactions on a square lattice, 
and arbitrary 1-local interactions, are QMA-complete. The first step of the proof develops a 
new variant of a reduction used in [12], in order to reduce such Hamiltonians to an XY model 
Hamiltonian on a spatially sparse graph. A key technical ingredient in the gadget construction 
is a careful analysis of the exactly solvable antiferromagnetic XY model on a cyclic chain. The 
second step is to apply variants of gadgets used by Oliveira and Terhal [22] to reduce the sparse 
graph to a lattice. Finally, to prove QMA-hardness of other interactions, we define new gadgets 
for reductions from the XY model which fit onto lattices. 

1.3 Organisation 

We begin, in Section 2, by describing the ideas from perturbation theory and the gadgets that 
we will use throughout the paper. Section 3 contains our results on interactions with restricted 
signs and StoqMA-completeness. Section 4 describes the QMA-hardness proofs for 2D lattices, 
with calculations regarding the XY model on a cyclic chain deferred to Section 5. We conclude 
with some open problems in Section 6. 
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2 Perturbation theory 


2.1 Definition of simulation 

We typically start with a target Hamiltonian Htarget acting on an Wdimensional Hilbert space 
of a form which is already known to be QMA-complete. The idea is to use a simulator 
Hamiltonian which acts on a larger Hilbert space Hsimj such that when restricted to its 
low energy subspace C]sf{Hsira) - that is, the subspace spanned by the N eigenvectors of Hgim 
with lowest eigenvalues - Hsim is approximately the same as the target Hamiltonian -Htarget- 
We formalise this idea of simulation, using the definitions presented in [8]. 

We require that there exists an isometry S ^ T-Lsim such that the image of £ is exactly 
C]\[{Hsira) and £^HgimS approximates Utarget in operator norm up to a small error e. This is 
enough to show that the first N eigenvalues of Hgim approximate the eigenvalues of H up to 
error e. 

However, this gives us no information about the eigenvectors of and how they relate to 
the eigenvectors of Htargetj as the isometry £ is unknown and may in general be very complicated. 
It will therefore be useful to have a simple isometry £ : % ^ T-isira that approximates £ in 
operator norm. This motivates the following definition of simulation as put forward by Bravyi 
and Hastings [8]: 

Definition 1. Let H be a Hamiltonian acting on a Hilbert space H of dimension N. A Hamil¬ 
tonian Hsim and an isometry (encoding) £ :H ^ Hgim are said to simulate H with error {r],e) 
if there eists an isometry £ : H ^ Hsim such that 

51. The image of £ coincides with the low-energy subspace CN{Hsim)- 

5 2 . \\ H -£^ Hsim £\\ ^ e- 

5 3 . \\£ — £\\ ^ rj . 

It follows from Weyl’s inequality [8] that, if {Hsim,£) simulates H with error (r/, e), the i’th 
smallest eigenvalues of Hsim and H differ by at most e for all I < i < A^. In addition, it is 
shown in [8] that simulation makes sense under composition as long as the lowest N eigenvalues 
are separated from the rest of the spectrum by a relatively large gap. 

Lemma 1 (Bravyi and Hastings [8]). Suppose {Hi,£i) simulates H with error {r]i,€i) and 
{H 2 ,£ 2 ) simulates Hi with error {r] 2 H 2 )- Tat Ai he the spectral gap separating the N smallest 
eigenvalues of Hi from the rest of the spectrum and suppose that Ai > 2 e 2 and ei,e 2 ^ ||-f^||- 
Then {H, T 2 T 1 ) simulates H with error (q, e), where q = qi q 2 + 0{e2Af^) and e = ei + 62 + 
0{e2Af^\\H\\). 

If Ai » ||H||, we have r/ 771 + 772 , e ~ ei + 62 . 

2.2 Perturbative gadgets 

The idea of a perturbative gadget is to have a simulator Hamiltonian consisting of two parts 
Hsim = AHq + V, for some large parameter A 3> 1, where ||H|| ^ A/2, so that the AHq term 
dominates. 

Let H- he the N dimensional ground space of Hq, and split the simulator space as Hsim = 
H- © H+. The encoding map £ : H ^ Hsim will always be chosen such that Im(T) = H_, and 
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will usually be obvious in context. Let P± be the projection operators onto 'H±. For any linear 
operator O acting on T-Lsim we define the following subscript notation: 

O— = P-OP- 0+ = P+OP+ 0-+ = P-OP+ 0+- = P+OP-. 

We assume that {Hq) _= 0 and that all eigenvalues of (i/o)++ are greater than 1. 

To find the effective Hamiltonian, we use the Schrieffer-Wolff transformation [7, 8]. This is 

a unitary operator on %siraj where S is an antihermitian operator satisfying S _= 0 = 

relates the subspaces and P-, such that {e^ _ \- = 0 = (e'^i7sime“‘^)+-. 

These equations in fact determine the operator S uniquely. Let £ = e~^£ so (SI) is satisfied by 
the definition of e~^, and for (S3) \\£ — £\\ = ||/ — e“‘^|| = 0(||5||) which it is possible to show 
is small for large A. 

Then define the effective Hamiltonian = {e^Hsime~^) , so 

= £^e^Hsime-^£ = £^H,s£ 


and so (S2) becomes 

IlLTtarget £^Hs\m£\\ — ||LL:arget ~ ||L7target Heff|| 

where Lftarget = TL^targetT^ is the logical encoding of the target Hamiltonian /^target in the 
simulator space. 

The aim is therefore to find a simulator Hamiltonian iLsim such that approximates 
Fftarget- Calculating LTeff exactly is very difficult, but we can express it as a Taylor series in 
powers of V. Truncating this Taylor series at first, second and third order gives the following 
three lemmas respectively. For a more detailed explanation of the Schrieffer-Wolf transformation 
and proofs of all of these results, see [8] . 

Lemma 2 (First Order). [8] Suppose one can choose Hq, V such that the non-zero eigenvalues 
of Hq are all greater than or equal to 1 and 


Htarget —|| ^ £/2. 


Suppose ||H|| ^ A. Then Hsim = AFfo + C simulates Htarget with error {rj,e), provided that 
A ^ H(e“^A^ -I- t/“^A). 


Lemma 3 (Second Order). [8] Suppose one can choose HQ,Vmain,y extra such that (FFo)+-i- tias 
all eigenvalues greater than or equal to 1, {Vextra)-\ _= {Vextra) _h = 0, {Vmain) _= 0, and 


\Htarget (Vextra) -T (Vmain) —h-^O (Vmain )-\—|| ^ 


Suppose 11 11 ? IlFeajtr'all ^ A. Then Hgim — AFFq -j- a ^ Vmain T Vextra simulates Htarget with 

error [t], e), provided that A ^ f2(e“^A® + r]~'^A‘^). 

Lemma 4 (Third Order). [8] Suppose one can choose Hq, Vmain, Vextra, Vextra such that (FFo)-i-+ 

has all eigenvalues greater than or equal to 1, (Vextraf^ _= (Vextra) _h = 0, (Vextra)^ _= 

(Vextraf —h — 0? (Vmain) -— 0; 


\Htarget (Vextra) T (Vmain) — l-^Q (Vr 


mainJ -\—h 


HQ^Vrr^ 


h-Ke/2 


and 



— (Vmain)—\-^Q (Vmain)-{ - 


Suppose ||Hmam||) ||Fea:ira||) IlFextrall ^ A. Then Hgim — AFFg -|- A ^ Vmain T A Vextra T Vextra 
simulates Htarget with error (r],e), provided that A ^ H(e“^A^^ + ?7“^A^). 


7 


Lemmas 1-4 reduce the problem of simulating desired effective interactions to the problem 
of finding suitable operators Lfoj 14iain, lextra, 14xtra- To prove hardness results we will always 
take e,rj = 0(1/poly(n)), making the approximation error small enough that it can be disre¬ 
garded. We do not attempt to optimise (or even calculate precisely) the size of the weights 
A = 0(poly(n)) that occur in our reductions, but they will always be polynomial in n and 
easily computable. 


2.3 Parallel use of gadgets 


It is often required to implement a poly(n) number of perturbative gadgets in parallel across 
different subsets of qubits without increasing the coefficients to super-polynomial size. Here we 
describe (following ideas of [22, 7]) how this can be achieved in a quite general setting. 

Consider a system of n qubits partitioned into k + 1 disjoint subsets of qubits, labelled Si 
for i = 1, 2,... /c. We will apply heavy interactions to the first k subsets, in order to simulate 
an effective Hamiltonian on the remaining subset of qubits Sk+i- Let Hq ' denote a heavy 

interaction which acts only on Si, be an interaction term that acts non trivially only on 

(i) 

Si U Sk+i, and P_^ the projection operator onto the ground space of Hq. 

Now consider the Hamiltonian terms for the whole system Hq = H^'^ and V = 

The projection operator onto the ground space of Hq is given by For a first order 

perturbation it is immediate that V _= so applying Lemma 2, we see that 

AHq + V simulates P--) as long as A = H(e“^||H|p). 


Since Hes acts only on H -, we need only consider how such an operator acts on the ground 
space of the overall Hamiltonian Hq, which has basis elements that can be expressed as a product 
state where each is in the ground space of Hq\ 

For second order perturbation terms we need to calculate Since acts 

nontrivially only on subsets i and A: -|- 1, only will be elevated out of the ground space by 
and so: 






hj 


i 


For third order perturbation terms satisfying hI*! : 
the fth subset of qubits in and out of the ground space H.-, 


- 0, each term can only move 
so 


F_+Ho-V++Ho-V+_ 






So in each case, the effective Hamiltonian simulated by H = ^ and Hq = is the 

sum of the interactions simulated by each of the and separately. 


2.4 The basic gadget 

We now describe a second-order perturbative gadget which will be used throughout this paper. 
Let H = aXX + I3YY -|- ^ZZ where q;-|-/ 3>0, a-|- 7>0 and /3 -|- 7 > 0. Consider a system 
of n qubits labelled by i G {1, 2,... n} along with two extra ancilla qubits a and b. The heavy 
interaction term Hq is just the interaction H applied to the ancilla qubits and a constant term 
to ensure that {Hq) _ = 0. This could then be multiplied by a further constant term to ensure 
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that (i?o)++ has all eigenvalues greater than 1, as required for Lemmas 2-4. Then 

Ho = l[Hab + {a + /3 + 7 )/] = (a + /3)|T+)(T+| + (a + 7)l^+)(^+l + (/? + 7)l^")(^”l ( 1 ) 

where are the maximally entangled states defined by = 

|oo)±|ii) 

V2 ■ 

In this diagonal form, we see that the conditions a-|-/3>0, a-|- 7>0 and /3 -|- 7 > 0 are 
exactly what is required for to be the unique ground state of H. Furthermore it is easy to 
calculate the inverse of H, which will be necessary to determine the second order perturbation: 


Ho^ = 


1 


;|T+)(T+| + 


1 




1 


|$~)(4>"|. 


a-|-/3 a- 1-7 /5 + 7 

To dehne V, partition the qubits {1, 2,... n} into two sets A and B; and let 

V — ^ ^ ^iHia T ^ ^ f^jHjb' 

i£A j^B 


Note that since |T )ab is maximally entangled, V _= 0. We also need to calculate 14_+, 

which can be done term by term: 

= aW|T-)(cI)-Ub + i/3y,|T-)(ch+Ufe - 7Zi|T-)(T+Ufe = -{Hib)-+ 

To calculate _, simply note that for any Hermitian operator V, _= This is now 

enough to calculate the second order perturbations. First 

^ = {H,b)- + H^HH,b)+- 

These terms just contribute to an overall energy shift in the Hamiltonian. The following terms 
simulate more interesting interactions: 

{H,,).+H^\Hjb)+- = (- 
{H,a)-+H^\HjaU- = 

iHib)-+H^\Hjb)+- = 

We define 


a 


/3 + 7 


XiX. - 


/3' 


-Y^Y, - 




^ a + 7 * ^ a+ ,5' 




\ab 


a 


/3 + 7 

o? 

/3 + 7 


/32 7^ 

XiX, + X—YiYj + -^ZiZj | 4 '-)( 4 '- \ab 


a + 7 

32 


a + j3' 
,,2 


XiXj + J^YiYj + T^ZiZ, ) |4'-)(T-Ub 


a + 7 


cx 13 


~ o? 7 ^ 

H = -—XX + yy + -^^zz 

p + 7 a + 7 a + B 


(2) 


and will use this notation throughout the paper. Lemma 3 tells us that for large A, the Hamil¬ 
tonian AHo + simulates, up to an overall energy shift, the following target Hamiltonian: 


H 


target 


E 

i£A,j£B 


XifijHij y ^ XiXjHij y ^ nifjijHij 


i,j&A 


i,jeB 


The gadget is illustrated in Figure 2. We see that a positive interaction is simulated between 
any two target qubits connected to opposite ancilla qubits and a negative interaction is simulated 
between any two target qubits connected to the same ancilla qubit. 
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Figure 2: Interaction graph for the basic gadget. In this figure, as throughout the paper, 
vertices represent qubits and edges represent interactions. Vertices which are filled in represent 
qubits between which we would like to produce effective interactions, while vertices not filled 
in represent ancilla qubits. Thick lines represent heavy interactions, while thin lines represent 
weak interactions. 


2.4.1 Antisymmetric case 


We will also use a similar gadget for antisymmetric interactions H. A general 2-local antisym¬ 
metric interaction has the form XZ — ZX up to normalisation. We will use the following basis, 
in which XZ — ZX is diagonal: 


IV’o) = 


| 00 ) + | 01 )-| 10 ) + | 11 ) 


Hi) = 


| 00 ) + | 01 ) + | 10 )-| 11 ) 


H 2 ) = 


-| 00 ) + | 01 ) + | 10 ) + | 11 ) 


IV’s) = 


| 00 )-| 01 ) + | 10 ) + | 11 ) 


We consider a Hamiltonian with the same interaction graph as in Figure 2. As before, the 
heavy interaction term Hq acts on a pair of ancilla qubits labelled a and 5, where 

Ho = liXaZk - ZaXi, + 21) = + H2){H + 2|V’3)(V’3|. 

And as before the other terms will be of the form 

V = XiiXaZi - ZaXi) + ^ ^ij{XbZj - ZbXj). 

i£A j£B 


Note that {XaZi — ZaXi) _= 0 = {X^Zi — Z^Xi)-- for all i, so V__ = 0. To work out the 

second order terms, we first calculate 


{XaZi - ZaXi)+_ = Zi\'ll;2){'lpo\ab “ Vj|V’l)(V'OIab 

(XbZi - ZbXi)+_ = Zi\‘ll;i){'ljjo\ab + Vi|V'2)(V’o|a6 

So the second order terms are 


{XaZi - ZaX,)_ + H^\XaZi - Z,W)+- 


2|V'o)(V'o|ab = {XbZi-ZbXi)_+H^\XbZj-ZbXj)+_ 


which just contribute an overall energy shift, and 

{XaZi - ZaXi)_+H^\XaZ, - = {XiX, + Z,Zj)\^o) {^Po\ab 

{XaZi - ZaXi)_ + H^\XbZ, - ZbXj)+_ = {Z,Xj - WZ,)|V’o)(^/>o|ab 

which allows us to simulate an XX -I- ZZ interaction between qubits connected to the same 
ancilla qubit, and a ZX — XZ interaction between qubits connected to opposite ancilla qubits. 
That is, by Lemma 2, the Hamiltonian Hgim = AHq -|- simulates the target Hamiltonian 


Hi 


target 


XiHj{XiZj 




z,Xj) 


^ ^ XiXj{XiXj -|- ZiZj) — ^ ^ iii^j{XiXj -\- ZiZj). 
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3 Positive weights 

We now have all the ingredients we need to prove Theorem 2. 


Lemma 5. Given a purely 2-local symmetric interaction H, which can be written in the form 
aXX -\- I 5 YY + ^ZZ, the problem {i7 }+-Hamiltonian is either: 

i) QMA-complete, if a ft > 0, a j > 0 and /? + 7 > 0 

ii) in StoqMA, otherwise. 


Proof. First we show that for case ii) the Hamiltonian H is stoquastic, in the correct choice 
of basis, and hence that the problem {Hj^-HAMILTONIAN is in StoqMA. Indeed, by a local 
change of basis on each qubit, we can relabel X,Y, and Z interchangeably. Therefore without 
loss of generality we can assume that a ^ /3 ^ 7 . We have 

0 0 a-f \ 

—7 a f 0 
a ft —7 0 

0 0 7 / 


H = aXX + fYY + 7 ZZ = 


/ 7 

0 
0 

\ a- f 


Looking at the matrix elements of H, the only off-diagonal elements that could be positive 
are those equal to a + ft, but since we are not in case i) and a ^ /3 ^ 7 , we must have that 
a -|- /3 ^ 0. So all off-diagonal elements are non-positive, which is the definition of a stoquastic 
Hamiltonian. 


Now we consider case i). We show the problem is QMA-complete, by showing we can use 
the basic gadget from Section 2.4 to simulate any Hamiltonian which is a sum of terms of the 
form 


~ cr 

H = -AA + ^ 


/3 + 7 


a + 7 


-YY + 


T 


Of. ft 


ZZ 


with both positive and negative weights. Observe that, asa-|-/3>0, a-|- 7>0 and /3 -|- 7 > 0, 
H is a sum of at least 2 Pauli terms. The Local Hamiltonian problem was shown in [12] to 
be QMA-complete for any interaction of this form. 

For each desired interaction between two of the n target qubits we wish to simulate, we 
add a pair of ancillary qubits, resulting in a total of at most n -t- 2n{n — l )/2 = physical 
qubits. Suppose for example that we wished to simulate an interaction between target qubits 
1 and 2. Then we introduce a pair of ancilla qubits 012 and 612 and apply a heavy 
interaction between them. Then if we wish to simulate a positive interaction we use a second 
order perturbation term V = + -^^ 6122 ; or for a negative interaction we use V = Hia^.^ + 

Hai 22 - As shown in Section 2.4, this will simulate the interactions H 12 and —H 12 respectively. 

These gadgets can then all be applied in parallel, as discussed in Section 2.3. □ 


Of those interactions contained in StoqMA, some of them are trivially contained in P. In 
particular, this holds when jOO) is the ground state of H. In this case the state jOOO ... 0) is the 
ground state for any Hamiltonian of the form Htot = Ylij with energy (OOjHjOO) Xij. 

Since and ]'!'=*=) are always eigenstates with energies as in eqn. (1), the state jOO) is a 
ground state when both and |<h“) have the smallest energy. That is, when 

a -|-7 = / 3 -|- 7 ^ min( 0 , a + ft) 


II 







which is equivalent to 


a = f3 and 7 ^ — |a|- 

Therefore for any interaction H which is, up to normalisation and relabelling X, Y, Z, of the 
form a(XX + YY) — ZZ for |q:| ^ 1, the problem {ff}'''-HAMlLTONlAN is contained in P. 

3.1 Antiferromagnetic TIM 

It remains to prove the StoqMA-completeness part of Theorem 2. Before we do so, we prove 
a related but somewhat simpler result about the Transverse Ising Model (TIM). This model 
consists of local interaction terms of the form ZZ and X\ 

H = atjZiZj + PkXk 

i,j k 

It may be assumed that /3fc is non-negative for all k, by conjugating all qubits for which /?*, < 0 by 
Z. This leaves the ZZ part of the Hamiltonian unchanged, and flips the sign of the affected /3fc. 
The problem of finding the ground-state energy of this Hamiltonian was shown to be StoqMA- 
complete when both positive and negative signs are allowed for , by Bravyi and Hastings [5] . 
However Bravyi [8] also showed that if aij ^ 0 for all i,j, then there is a polynomial-time 
probabilistic classical algorithm to find the ground-state energy (and even to approximate the 
partition function), so the problem {— ZZ, A}'*'-Hamiltonian is in BPP. 

We therefore focus on the antiferromagnetic case, where aij ^ 0 for all i,j. 

Theorem 5. The problem {ZZ, A} -Hamiltonian is StoqMA-complete. 

Proof. Since the full TIM Hamiltonian has been shown to be StoqMA-complete it will suffice 
to show that we can simulate this Hamiltonian using only positive aij. We do this by encoding 
logical qubits in the two dimensional ground space of ZZ, which is span{|01), |10)}. We label 
logical states and operators with a L superscript, identifying |0'^) = IT”) and |1^) = I'b'*’), 
and associate each logical qubit i with physical qubits ia and %■ We will use operators of the 
following form to simulate and Z^. 


Ho — ^ 5 ^ |00)(00|j -I- 111) (11 |i 

i=l i=l 

14xtra = 

i,j:\ij>0 i,j:\ij<0 

f^nain = hi i^ia + ^ib ) 

i 

Calculating the first order perturbations gives: 

(Z,,Z,J__ = (|T+)(T-|i + |T-)(T+|,) 0 (|^&+)('b-|, + |^-)('b+|,) = XtX^ 

(Z,,Z,J__ = (|'b+)(T-|, + |'b-)(T+|,) 0 (-|T+)(T-|, - |'b-)(T+|,) = -XtX^. 

This is how it is possible to simulate both XX and —XX. Note that (A*^)_= 0 = (Aj,^) _ 

for all i, so I4iain has no first order contribution. To calculate the second order perturbations 
for Hmain, we first show 


(A,, + XO+_ = (|00) + |11))((01| + (10|), = 2|cb+)(T+|, 
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so 

(X,, + = 4|vl/+)(vl/+|, = 2{I - Zl) 

Therefore, by Lemma 3, the Hamiltonian -Hsim = + A^/^ Vmain + 14xtra simulates the target 

Hamiltonian: 

-^target — ^ ^ ^ij^iZj T E — Zi). 

i,j i 

Relabelling X and Z gives the general TIM Hamiltonian, which was shown to be StoqMA- 
complete in [5]. □ 


3.2 The case H = XX — YY + 'yZZ for 7 > 1 

We finally show that the remaining cases on the boundary of the StoqMA region are also 
StoqMA-complete, also using a reduction from TIM. The proof is similar to the previous section. 

In the case where a = 1 = —/? and 7 > 1, then both |4'“) and I'L'*') are ground states 
for the interaction H. Therefore we can take a system of 2n qubits labelled by ia and if, for 
i = 1,2,.. .n, and let 

n n 

Ho = Y,\ + 7/) = + (7 - 

i=l i=l 


Hextra — ^ijHia,jb 

i,j.Xij>0 i,j-.\ij<0 


4^ain — 


aijb / 




The encoding S encodes a logical qubit into each pair of qubits, with jO-^) = |T ) and |1^) = 

|T+). 

Calculating the perturbations to first order gives, similarly to before. 




= 7(l'i'+)('i^“k + ® - l^-)('k+|,) = -7XtX^ 

Therefore (I4xtra) _= (Hmain) _= 0. To calculate the second order 

contribution of fmainj consider just one term, V = 

H+_ = 2(|$+)|$-) - |$-)|4>+))(T-|(T+| + 2(|$+)|$+) - |$-)|$-))(T+|(^+| 


Since |<I>'*’)|‘h ) and |<h both have energy 27 with respect to the Hamiltonian Hq; and 

|(J)+)|(|)+) and |<I>“)|<h“) have energies 27 + 2 and 27 — 2 respectively, we get: 


H_+ifo-V+_ = 4|4/-)|T+)(T-|(T+| + ) 4|4/+)|T+)(T+|(T- 


L-,L\inL-,L\ , ^7 nL.L\/iLiL| 


= -|0^W)(0^W| + „ ^ 

7 7^ — 1 




7(7^ -1) 


{Z^ - (272 - 1)1) ® {Z^ - I) 
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Then by Lemma 3, the total Hamiltonian ffsim = Hq + Al/gxtra + A^/^T main simulates the 
Hamiltonian 

- A, 1) (Z - (27= - 1)'). (Z - 

and so we can implement any Hamiltonian which consists of interactions XX and ^( 7 ) = —{Z— 
( 27 ^ — 1)I){Z — I). Simulating a TIM Hamiltonian requires one further round of perturbation 
theory to simulate Z interactions. We add an extra pair of ancilla qubits with a heavy F(^) 
interaction applied between them to project them into the ground state |11). Then we apply 
F{^) between one of these ancilla qubits and another qubit, to simulate a Z interaction to this 
qubit. 

This completes the proof of Theorem 2. 


4 Restricted interaction graphs 

We now move on to proving QMA-completeness results for interactions with restricted geome¬ 
tries. First, we show that the XY model is QMA-complete even when the interactions are 
restricted to the edges of a triangular lattice, and have positive weights. To prove this we first 
use similar methods to [ 12 ] to show that the model is QMA-complete on a spatially sparse 
graph, as defined by Oliveira and Terhal in [22]. Then we generalise the fork and crossing 
gadgets used in [22] to reduce this to a 2D triangular lattice. 

4.1 Spatial sparsity 

In [12], it is observed that the results of [22] and [3] can be straightforwardly combined to show 
that a Hamiltonian of XX, ZZ, X and Z terms is QMA-complete even when the interaction 
terms are restricted to the edges of a 2D square lattice. We will simulate such a Hamilto¬ 
nian using a variant of the gadget used in [12] to show that the XY model is QMA-complete. 
Unlike [12], our gadget will fit onto a spatially sparse graph. 

Definition 2 (Oliveira and Terhal [22]). A spatially sparse interaction graph G is defined as a 
graph in which (i) every vertex participates in 0 ( 1 ) edges, (ii) there is a straight-line drawing 
in the plane such that every edge overlaps with 0 ( 1 ) other edges and the length of every edge is 
0 ( 1 ). 

Let H = XX + YY. In the construction of [12], three physical qubits are used to encode 
two logical qubits in the ground space of H 12 + Ft 23 -I-LI 13 , which is the asymmetric subspace of 
three qubits. It is then shown that XI and ZI can be simulated by applying further interactions 
within this subspace, and that applying interactions across triples of physical qubits can result 
in logical interactions of the form: 


X,Xj{Xi,Xp + Yi,Yj,) 


Z,Zj{Xi,Xy + YeYj,) 
hIj{Xi,Xj, + Y,,Yj,) 

where we have labelled the first and second qubit in each pair of logical qubits as i and i' 
respectively. We can use the third type of interaction term to project all of the second qubits 
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in each pair into the ground state of a certain Hamiltonian Hget- Having done this we can use 
the first two terms to create a Hamiltonian of the form: 

+ PkZk) + + 6ijZiZj){n\Xi,Xj, + Yi,Yj,\n) 

k i<j 

where |n) is the ground state of Hget- 

Therefore it is possible to implement any Hamiltonian consisting of terms {XX, ZZ, X, Z}, 
as long as we can find a Hamiltonian Hget such that: 

1. The ground state |n) is unique; 

2. The spin correlation functions {Q\XiiXj/ + YiiYjr\Q) can be calculated in poly(n) time on 

a classical computer, so that the coefficients can be suitably adjusted; 

3. In order that are 0(poly(n)), it is also necessary that the correlation functions are 

n(l/poly(re)); 

4. Similarly we need the spectral gap of to be n(l/poly(n)). 

In the original proof the XY model on the complete graph was used [12], but this does 
not satisfy the spatial sparsity constraints. To meet this constraint, we instead consider a ID 
model. In this context simultaneously satisfying conditions (3) and (4) is challenging. We 
require D(l/poly(n)) size correlation functions for qubits at large distance apart on the chain. 
But, as it is known that for any system of this form with a constant spectral gap the correlation 
functions will decay exponentially [15], we require that the spectral gap —?• 0. However, by (4) 
we also require that the gap does not go to zero too quickly. 

Fortunately, however, the XY model on a cyclic chain satisfies all of these constraints: 

Lemma 6. Fix N even but not a multiple of 4, let H = '^^^^{XiXi^i + YiYi^i)+XiX]y + 
YiYat. Then H has a nondegenerate ground state jD) and spectral gap D(l/Y). Further, for 
any pair i, j such that \i — j\ = n and n = o(Y'^/’^), {Q\XiXj + YfYjl^) = D(n" There is 

an efficient classical algorithm to compute the spectral gap and all the correlation functions. 

The proof of Lemma 6 is deferred to Section 5. Most of the claims in the lemma are well- 
known [18, 19]; however, we were unable to find in the literature a rigorous lower bound on the 
correlation functions for a finite-length chain which is as tight as we need. 

To simulate a Hamiltonian H which consists of n qubits on a square lattice, we first add 
^/n ancilla qubits, which are not involved in any interactions, along one edge of the lattice (see 
Figure 3). Then, for each of the qubits in this extended lattice, we use three physical qubits 
to encode two logical qubits and apply the Hamiltonian of the XY model on a cyclic chain to 
the second qubit in each pair with the path of the cycle arranged as in Figure 3. Because the 
ancilla qubits did not interact with any others in the original lattice, the longest distance along 
this path for which we need to calculate the correlation functions (in order to simulate XX and 
ZZ interactions) is 2n^/^. By Lemma 6, the correlation functions are D(n“^/^). 

We have now obtained a spatially sparse Hamiltonian which only uses XY interactions and 
simulates H. However, the interactions no longer take place on a square lattice, because each 
site on the lattice contains 3 physical qubits in a triangle, and there are complicated interactions 
between neighbouring triangles. To go back to a lattice geometry we need some further gadgetry. 
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Figure 3: Cyclic XY model Hamiltonian overlaid on a 2D square lattice, with ^/n ancilla qubits 
(coloured white). 

4.2 Mediator qubit pair gadgets 

In [22], Oliveira and Terhal develop “subdivision”, “fork” and “crossing” gadgets involving a 
mediator qubit. This is an ancilla qubit a, which is projected into the state |0) by a heavy I — Z 
interaction, such that when 2-local terms are applied, such as F = Hia + H 2 a, second order 
perturbation theory gives an effective interaction between qubits 1 and 2, even though there is 
no direct interaction between them. This allows one to simulate a general Hamiltonian on a 
spatially sparse interaction graph using a Hamiltonian on a planar 2D graph. 

However, here we only have access to 2-local terms. The direct analogue of these gadgets is 
to take a pair of qubits, a and b and project these into the ground state of a 2-local interaction 
Hq = Hab- The second order perturbation theory analysis of this situation is exactly the basic 
gadget described in Section 2.4; and so the interactions generated will be of the form H (see 
eqn. (2)). We focus on the XY interaction, where H = XX + YY and so H = H (this is similar 
for the Heisenberg interaction H = XX + YY + ZZ, so all of these gadgets would work in 
this case too). For the fork and crossing gadgets, where there are three or more extra qubits, 
unwanted interaction terms are generated that need to be cancelled out by first order terms. 

4.2.1 Subdivision gadgets 

The simplest gadgets are the subdivision gadgets, which are in a sense the most basic gadgets 
we can make using the XY interaction. The interaction graphs for the gadgets are shown in 
Figure 4. There are two additional qubits, labelled 1 and 2, as well as the mediator qubits, 

b O 

•-0—0-• •-O-• 

1 a b 2 1 a 2 

(a) Positive subdivision gadget (b) Negative subdivision gadget 

Figure 4; Subdivision gadgets 
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(a) Fork gadget (b) Crossing gadget 


Figure 5: Fork and crossing gadgets. In each case the left-hand interaction pattern is simulated 
using the right-hand gadget. 


which are either connected to the same or opposite ancilla qubits. 

When connected to opposite ancilla qubits, i.e. V = Hia + H 2 b, the analysis of Section 2.4 
shows that AHq -|- A}/‘^V simulates a X 1 X 2 + Y 1 Y 2 interaction. And when connected to the 
same ancilla qubit, i.e. V = Hia + H 2 a, the analysis of Section 2.4 shows that AHq -|- A^/'^V 
simulates a —X 1 X 2 — Y 1 Y 2 interaction. 

The positive subdivision gadget can be applied in series, that is it can be used to simulate any 
of the interactions Hia, Hat, and Hb 2 , using a new heavy weight A' ^ A to give a Hamiltonian 
on 6 qubits in a line, with an effective XX + YY interaction between the hrst and last qubits. 
This can be repeated at most a constant number of times giving a Hamiltonian on 2k qubits 
on a line, where k = 0(1), which simulates an XX + YY interaction between the qubits at the 
endpoints. 

Similarly for the negative gadget, both Hia and H 2 a can be simulated by the positive 
subdivision gadget, and this process can be repeated 0(1) times. Thus for any k = 0(1), there 
exists a Hamiltonian on 2k qubits where 2fc — 1 qubits are on a line, and one extra qubit is 
connected to a qubit at an odd internal position on this line, which simulates a —XX — YY 
interaction between the qubits at the endpoints of the line. 

4.2.2 Fork gadget 

The fork gadget is used to reduce the degree of a vertex in an interaction graph. For a qubit 
interacting with two other qubits, we can simulate this using a gadget which only has one 
incoming edge to this qubit. The interaction graph is given in Figure 5a, and has I4iain = 
Hia + H 2 a + H^b- This will also generate an unwanted —Hu, so it is necessary to also include 
the term I4xtra = Hu, which is represented on the interaction graph by a dotted line. Then, 
by Lemma 3, the Hamiltonian Hgim = AHq + A^/^I4iain + lextra simulates Hu + H 23 . 

4.2.3 Crossing gadget 

The crossing gadget is used to simulate a Hamiltonian in which two edges in the interaction 
graph cross, using a gadget with no crossings, as can be seen in Figure 5b. The main second 
order terms are Hmain = Hu + H 2 a + H^b + H^b: and the required correction terms are I4xtra = 
Hu + H 34 — Hu — H 24 . Then, by Lemma 3, the Hamiltonian Hgim = AHq + A^/^Fmain + I4xtra 
simulates Hu + H 23 . 
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4.3 XY model on a 2D planar graph of degree 3 

We can now follow the prescription of Oliveira and Terhal [22], using the analogous mediator 
qubit gadgets described in Section 4.2, to simulate a spatially sparse Hamiltonian using a 
triangular lattice. This follows [22] closely so we only sketch the steps involved. 

First use the subdivision gadget on every single edge in parallel in order to isolate each 
vertex of high degree. Then, for each vertex of degree d = 0(1), pair up adjacent incoming 
edges and use the fork gadget in parallel on each pair of incoming edges, to reduce the degree 
of the vertex to [d/2]. Then repeat this process O(logd) times in series until the vertex has 
degree at most 3. Each stage of this process can be done in parallel across all vertices of high 
degree, such that only a constant number of reductive steps are required to end up with an 
interaction graph of degree 3. 

Then we isolate each crossing by applying the subdivision gadget up to O(logm) times on 
each edge in series, where m = 0(1) is the maximum number of crossings across any edge in 
the graph. Finally we apply the crossing gadget to every crossing in parallel to be left with 
a Hamiltonian with a planar interaction graph G of degree at most 3, such that there is a 
straight-line drawing in the plane where each edge is of length 0(1) and the angle between 
adjacent incoming edges at a vertex is 0(1). 

4.4 2D triangular lattice 

We finally simulate the Hamiltonian of the previous section using a Hamiltonian on a 2D 
triangular lattice with only positive weights. 

To do this we take the interaction graph G from the previous section, along with its em¬ 
bedding in the 2D plane such that the graph is planar with maximum degree 3, each edge has 
length 0(1), and the angle between adjacent edges is 0(1). Then we place a fine triangular grid 
on top of this graph and move each vertex of G onto the nearest lattice point, then deform edges 
of G onto the nearest path through the lattice. Since the angle between any two edges is 0(1), 
if the spacing in the lattice is sufficiently small (but still 0(1)), these paths will only intersect 
in a small 0(1) region around each original vertex. The paths can then easily be rerouted in 
this area. 

As discussed in Section 4.2, the subdivision gadget can be used in parallel to simulate +H 
(resp. —H) when there are an even (resp. odd) number of mediator qubits on the path in 
between the target qubits. On a triangular lattice it is easy to include one extra qubit on the 
path, for example by taking two edges of a triangle rather than one, and so it is possible to 
choose the appropriate parity, and simulate either +H or —H as desired. This completes the 
proof that the antiferromagnetic XY interaction on a triangular lattice is QMA-complete. 

Note that this construction will work for essentially any non-bipartite lattice. A graph is 
non-bipartite if and only if it contains an odd cycle. If the spacing of the lattice is small enough 
such that there is such an odd cycle between each connected pair of original vertices, and this 
odd cycle is connected to the rest of the graph at two points, then two paths that pass through 
these two points (but which go around this cycle in opposite directions) will have lengths of 
opposite parities, and thus we can make either +H or —H as desired. 

4.5 Other interactions on a triangular lattice 

We now return to considering more general interactions of the form H = aXX+l3YY-\-^ZZ. In 
order to prove that these interactions are QMA-complete on a triangular lattice, we will simulate 
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(a) Gadget on a triangular lattice (b) Hamiltonian simulated by this gadget 

Figure 6: Triangular lattice gadget for simulating XX + YY 


a target Hamiltonian of XY interactions on a triangular lattice using only local gadgets. We 
require that when each edge of the target interaction graph is replaced by the local gadget, the 
resulting graph (the interaction graph for the simulator Hamiltonian) is still a triangular lattice. 
Therefore each local gadget must fit on a triangular lattice, and in order that separate gadgets 
do not interfere, they must be able to tessellate without overlapping. Figure 7 shows that the 
triangular gadget described below does indeed tessellate in this way. 

The idea behind the gadget is to simulate two different interactions and along 
two different paths between two qubits, and then to take a linear combination of these to access 
a different interaction known to be QMA-complete, such as the XY interaction. Figure 6 (b) 
shows such a gadget that will simulate two different interactions between qubits 1 and 6 (where 
H is as in (2) and H' is defined below). However, both qubits 1 and 6 have two incoming edges, 
and so this gadget will not tessellate on the lattice. We therefore need to use an equivalent of 
the fork gadget so that qubits 1 and 6 only have one incoming edge each. 

The interaction graph of the triangular lattice gadget is shown in Figure 6 (a), and it will 
simulate the Hamiltonian shown in Figure 6 (b). The red and blue sections of the graph are 
generalizations of the Fork gadget for interactions of the form H = aXX+j3YY+^ZZ. Each of 
them consists of two parallel applications of the basic gadget. Consider the fork gadget coloured 
blue in Figure 6 ; the mediator qubit pair connected to qubit 1 generates interactions of the form 
+Hi 2 and +H 14 , but also an unwanted interaction term —H 2 A between qubits 2 and 4. This is 
cancelled out by simulating a +H 24 interaction with the other pair of mediator qubits between 
qubits 2 and 4. If both of the physical H interactions acting on qubit 2 in this fork gadget have 
weight A, and both interactions acting on qubit 4 have weight /x, then the overall effect of the 
blue fork gadget is to simulate AH 12 + /XH 14 . 

Similarly the part of the gadget coloured red can be viewed as a kind of fork gadget. One me¬ 
diator qubit pair simulates and - 1-7735 interactions, while the other pair of mediator 

qubits simulates — 7/35 resulting in an overall effective Hamiltonian of — H^q. 

The gadget between qubits 4 and 5 is just the basic gadget and so simulates H 45 , but 
the gadget between qubits 2 and 3 is slightly different and needs to be studied separately. 
Labelling the qubits in between 2 and 3 as a,b,c, d, the heavy Hamiltonian for this gadget is 
Hq = Hab + Hcd + 2(a + j3 + 7)7 which has non-degenerate ground state The 

other terms are V = 772 a + Hbc + Hds, and it is straightforward to show that 

F__ = 0, and V-+Hq^V+- oc 7 

To calculate the third order contribution, it will help to write out V as nine terms, 

H = aY2Ya + mYa + lZ2Za + aXbX^ + ^YbYc + + aX^X^ + (dY^Ys + jZ^Z^, 

and work in the basis I'k^), for the qubit pairs (a, b) and (c, d). Then each term of V maps 
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basis elements to basis elements, and Hq is diagonal. 




a 


7 X 2 X 3 + 


/3' 


7X2X3 + 


7 


7X2X3 


(/3 + 7)2 (a+ 7)2 (q! + /3)2 

Therefore, by Lemma 4, the Hamiltonian AL^o + simulates the interaction H'\ 


H' = 


(/? + 7) 


2 (a+ 7)2 (a + Z?)^ 


zz =: a A A + p 


Note that the coefficients of H' satisfy a' + /3' > 0, /3' + 7 ' > 0 and a' + 7 ' > 0, so LT' can 
be used as the heavy interaction term Hq in a gadget similar to that of Section 2.4. 

Therefore, overall the gadget simulates the Hamiltonian shown in Figure 6 (b). This Hamil¬ 
tonian is then used to simulate an interaction (via qubits 2 and 3) and an interaction 

(via qubits 4 and 5), where 


a 




7 


— 7 XX -h —-YY + ^, ZZ 

/?' -|- 7' a + 7 a + fi 


iL( 2 ) = 


a 


/3 + 7 


-XX + 


/32 


a -I- 7 


zYY + 


7 


a. Y (3 


.ZZ. 


We now show that having access to interactions of the form and H^‘^\ with arbitrary 

positive weights, is sufficient to produce an AX -|- YY interaction. 


4.5.1 The case aXX + I3YY 

First we consider the case where H has Pauli rank 2, where 7 = 0, and so a, /3 > 0. In this case 
it is particularly easy to explicitly find S, a', (3, [3' in terms of a and (3 to show that 

= + and if® = iZ-y X + Wr, 

(3^ a° [3^ 

In this case we can directly simulate the XX -|- YY interaction with Y by taking 

a® — 

a^/3'^(a2 — /32) ^ a®/3®(a2 — /32) 

The exact values that A and p take are not too significant; but it is important to note that they 
are both positive (given that a, /3 > 0) and are easily computable. 

4.5.2 The case aXX -|- ^YY Y 7 XZ 

For more general interactions FI, where all of a, /3 ,7 are non-zero, the simulated interactions 
and FI( 2 ) ^re guaranteed to be different unless H is the Heisenberg interaction XX -|-Xy -|- ZZ. 
In order to normalise —let A = 1 — p, and define a(p),/3(p), 7 (p) to be the 
coefficients of 

H{n) = 

For p G [0,1], this is the one parameter family of interactions that can be simulated using the 
gadget. So H{0) = which has all coefficients negative and H{1) = which has all 
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coefficients positive. Therefore there exists fia, f^'y £ (0) 1) such that a{fia) = 0,= 
Oj 'y{fJ"y) = 0 - In particular: 

Q.(l) ^(1) /y(l) 

= a(i) + «( 2 ) ’ = /?(!)+/ 3 ( 2 )’ = 7 ( 1 ) + 7 ( 2 ) 

Calculating the coefficients q;(^),/ 3(/i), 7 (^) exactly in terms of a,/ 3,7 and // yields very 
messy expressions, but the following lemma gives a very useful relation between a, (3 ,7 and 
h'aj IJ'fii fJ‘'y 

Lemma 7. If'j = a, then = fj-a- If 1 > oi and 7^/3, then > ^a- 


Proof. The aim is to show that — Ha ^ Q- First simply substitute in the expressions for n.y 
and ^a- 

7(1) Q,(i) 0:07^) — q-(i)7(2) 

7(11-1-7(2) Q,(l) -g q,(2) (1^(1) _|_ ^(2)^^q,(1) _|_ q,(2)^ 

Since every term in the denominator is positive, it will suffice to consider just the numerator. 

(/?' + 7 ')(5 + ^) - (^ + 7 )(o' + j3') 


q,( 2 )^( 1 ) _ Q,(l).y( 2 ) — 


5272 


(^' + 7')(/3 + 7)(«' + ^')(« + 


Again, the factor outside the square brackets is strictly positive, so we can just consider the 
expression inside the square brackets. The following relation will be useful; 


fl'a — [3a' = 


a2^2 


/3 


a 


a^f3^{a + /3 + 7)(/3 - a) 


(a + 7)(/3 + 7 ) ya + 7 l3 + jJ (a + 7)2(^-j- 7)2 

The expression in the square brackets {/3'a — [3a') + {'^'a — 70 ') + ( 7 'a — 70 ') is then just 
a + /3 + 7 


(a + 7 ) 2(/3 + 7 ) 2 (a + ,5)2 


[a^^‘^{a + /3)2(/3 - a) + l3^-f‘^{l3 + 7 ) 2(7 - /3) + a^j^{a + 7 ) 2(7 - a)] 


Note that setting a = 7 at this point would give zero, implying = /Xq,, thereby proving the 
first part of the Lemma. Now considering the case where 7^/5 and 7 > a, we can use the 
inequality [3'^^‘^{[3+ j)‘^{'y — /3) ^ [3“^a^{[3 + a)"^— [3) to replace the second term and show that 
the previous line is greater than or equal to 


a + /3 + 7 


{a + 7 ) 2(/3 + 7 ) 2 (a + /3)2 
which is strictly positive since 'j > a 


[q;2/32(q- -I- /3)2(7 - a) + o?'^‘^{a + 7)^(7 - a)] 


□ 


There are then two cases we need to consider: 


(i) a, /3, 7 have a unique maximum, say j > a, [3. 

Then by Lemma 7, > iXa,hi3^ so a{ix^) and [3{ix.y) are both positive and we can 

simulate II{fx.y) = a{ix^)XX + [3{ij,^)YY which is QMA-complete on the triangular lattice 
by Section 4.5.1. 

(ii) a < [3 = J. 

Then by Lemma 7, fXa < fx/s = ix^, so we can simulate II{^^) = a{fx.y)XX. So for e > 0 
small enough, fx = ix.y + s satisfies a{fx) > /3(/i),7(/x) and a(/x),/3(^), 7 (/x) all positive. So 
H (/x) is of the form we have just shown to be QMA-complete on the triangular lattice in 

(i). 
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Figure 7: Tessellations of the square and triangular lattice gadgets. 


This completes the proof of Theorem 4. 

Theorem 4 (restated). Let H = aXX + jSYY + ^ZZ be a 2-qubit interaction such that 
a + /3>0, a + 7>0,/3 + 7>0 and H is not proportional to XX + YY + ZZ. Then 
Hamiltonian is QMA-complete, even if the interactions are restricted to the edges of a 2D 
triangular lattice. 


4.6 2D square lattice 

A 2D square lattice is a bipartite graph, and for any bipartite graph we can locally change basis 
for every qubit in one side of the partition by conjugating by Z. Restricted to positive weights, 
this will take any QMA-complete interaction to one contained in StoqMA, by effectively taking 
a, /3 —>■ —a, —/3. Since we do not expect QMA=StoqMA, we will instead study the case where 
both positive and negative weights are allowed. 

To show that the XY interaction is QMA-complete on a 2D square lattice, we proceed as 
with the triangular lattice. We place a fine square grid over the 2D planar interaction graph 
of Section 4.3, move vertices to the nearest lattice points and deform edges to paths on the 
lattice (re-routing in a small region around each vertex if necessary to avoid collisions). We 
then use the subdivision gadget (Section 4.2.1) to implement effective interactions across the 
ends of the paths. However, on a bipartite graph such as the square lattice, the parity of the 
length of any path between two points is fixed. Therefore, if there are an odd number of qubits 
on the path but we wish to simulate a positive interaction (or vice versa), we need one of the 
interactions along this path to be negative. This is then enough to show that the XY interaction 
is QMA-complete on a 2D square lattice with positive and negative weights. 

For a more general symmetric interaction H = aXX-\-l3YYY^ZZ, given that we are allowed 
both positive and negative weights and can conjugate one side of the partition of the graph by 
X, Y or Z, we can assume without loss of generality that a, /3,7 are all non-negative. Figure 8 
shows a gadget that fits on the 2D square lattice and simulates the interaction 
in exactly the same way as the gadget on the triangular lattice, by first simulating the six qubit 
Hamiltonian shown in Figure 6(b). Again the blue and red gadgets act as fork gadgets, but 
in order for the unwanted H 23 and interactions to be correctly cancelled out, the sign of 
the interactions marked as dashed lines in Figure 8 must be negative. The other dashed line in 
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Figure 8: Square lattice gadget for simulating XX + YV. Dashed lines represent negative 
interactions. 


Figure 8 must also be negative so that the gadget between qubits 2 and 3 correctly simulates 
+^ 23 - It I® shown in Figure 7 that this gadget will tessellate on a square lattice without 
overlapping. 

Therefore, using the prescription set out for the triangular lattice in the previous section, 
we can use to generate XX + YY. We have thus shown that the interaction 

aXX + PYY + 'yZZ is QMA-complete on a square lattice with positive and negative weights, 
as long as two of a ,/?,7 are non-zero, and they are not all equal. 

4.6.1 Interactions with non-trivial 1-local part 

We will finally consider the more general case of 2-qubit interactions with a non-trivial 1- 
local part. It is sufficient to assume that 77^“^* is symmetric or antisymmetric under interchange 
of the two qubits [12]. For a symmetric Hamiltonian we can further assume that 

= aXX j3YY + jZZ + AI + IA = H + AI + IA 

where H is a general 1-local term. An antisymmetric interaction has, up to normalization, the 
normal form XZ — ZX + AI — I A. 

In both symmetric and antisymmetric cases, it is shown in [12] that the Hamiltonian 

u ufull I zrfull _ u zj \ u zj 

^0 — J^ab - ^cb + ^cd - ^ad — ^ab “ J^cb + ^cd “ ^ad 

has a unique ground state jH), and that in this state jH) the reduced density matrix for each of 
the qubits a, b, c,d is I / 2 . 

Therefore if we project into the ground state of the Hamiltonian Hq, and apply an extra 
interaction V = for example, then the hrst order perturbation term will be 

v.. = {n\-Hlf\n) = -Ae. 

So by Lemma 2 , we can simulate any multiple of —A^, which we could use to cancel out the 
1 -local part of the interaction acting on e. 

This gadget can be implemented on a square lattice as shown in Figure 9, to simulate 
arbitrary multiples of the 1-local term A on each of the qubits connected to the central square. 

These gadgets can be applied in parallel over a square lattice divided into 3x3 blocks, with 
additional physical interactions along the outside of each block, to simulate a Hamiltonian 
of the form depicted in Figure 10. All of the physical interactions shown on this interaction 
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Figure 9: Square lattice gadget for removing 1-local terms 


graph are of the form but the 1-local terms simulated using these gadgets are chosen such 

that they exactly cancel out the 1-local part of on all the white qubits. 

We label each black qubit with a label i, and label the pair of white qubits between i and 
j as aij and aji, where aij is the qubit nearest to i. We will construct a Hamiltonian acting on 
the interaction graph of Figure 10 that simulates a general Hamiltonian of the form ^ XijHij, 
where interactions take place between adjacent vertices on a sqnare lattice, which was shown 
to be QMA-complete in the previous section. 


On each pair of white qnbits we apply a heavy 2-local H interaction, so that the overall 
heavily weighted term is 

^ J 


Hn = 


^ |A, 


rH, 


ijl 


^ij ^ji 


where E is the set of edges of the sqnare lattice and J is a large weight to be defined. 

We choose the weights of the other interactions so that the 1-local terms all cancel out on 
the black qnbits, and then we can simulate H interactions as nsual by applying the basic gadget 
to the resulting H interactions. In order to do this we pick all the other weights (those between 
one black and one white qnbit) to be of the same magnitude \/J, where J = 0(poly(n)) is the 
max size of the weights in the target Hamiltonian, and let the sign of the interaction between 
qnbits i and aij be given by mj = ±1, so overall we have 


y — fJ'ijiHiaij + Hj) -|- + Aj). 


Therefore, ensnring that the 1-local terms cancel out on the ith black qnbit is equivalent to 
choosing Hij snch that Ylj f^ij = 0. In order to simulate the target Hamiltonian Htarget = 
X; XijHij, we need mj^ji = sgn(Ajj). 



Figure 10: Hamiltonian simulated by parallel use of gadget 
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(a) Basic gadget (b) 6 qubit gadget 


Figure 11: Gadgets for simulating XX + ZZ using XZ — ZX interactions. The 6 qubit gadget 
simulates the basic gadget to implement an effective XX + ZZ interaction across the left- and 
rightmost qubits. 

Suppose the black qubits along the top row are labelled 1,2,... from left to right. Pick 
(arbitrarily) ^12 = -|-1, and then choose inductively all the /ijj values along the top row according 
to the following rules: 

l^i+i,i — sgn(Aj^j-|_i and 

Then do the same for every other row as well, and repeat this procedure for each column. Then 
only the black qubits on the edge of the lattice still have 1-local terms, and these can be removed 
simply by adding extra square gadgets to the outside of the lattice. 

Therefore in the symmetric case, where H = aXX -\- j3YY Y'jZZ, the overall Hamiltonian 
AHo -|- will simulate the effective Hamiltonian 

Heff =-H.+Hq-V+_ = XijHij, 

(*d)6E 

where H is defined as in eqn. (2). Thus we can simulate any Hamiltonian consisting of H terms 
on a square lattice. In the antisymmetric case, H = XZ — ZX, we can choose signs in exactly 
the same way so that all 1-local terms are cancelled out and we simulate the target Hamiltonian 

^target — ^ ^ Xij {^X^Zj Z^Xj') . 

It remains to show that this Hamiltonian is QMA-complete, even when restricted to a 2D 
square lattice. We achieve this by simulating XX + ZZ terms on the lattice. When two 
qubits are connected to the same qubit of a mediator pair using XZ — ZX interactions, as in 
Figure 11a, we showed in Section 2.4.1 that an effective XX + ZZ is simulated. However, this 
gadget does not tessellate on the lattice, as there is not enough space to fit the two mediator 
qubits between physical qubits. In Section 2.4.1, we also showed that two qubits connected to 
opposite ends of a mediator qubit pair can simulate an effective XZ — ZX interaction, so this 
is a kind of subdivision gadget. If one of the edges of the gadget in Figure 11a is simulated 
using this subdivision gadget, we obtain the 6 qubit gadget in Figure 11b. This 6 qubit gadget 
will tessellate on a 2D square lattice, so we can use it to simulate the XY interaction (after 
relabelling Z to Y) on a 2D square lattice. 

We have finally completed the proof of Theorem 3. 

Theorem 3 (restated). Let H he a 2-local qubit interaction with Pauli rank at least 2, sueh 
that the 2-local part of H is not proportional to XX -|- YY + ZZ. Then {H}-Hamiltonian is 
QMA-complete, even when the interactions are restricted to a 2D square lattice. 

5 The XY model on a cyclic chain 

This section contains the deferred proof of Lemma 6. Recall that in this lemma we consider 
the XY model on a cyclic quantum spin chain of N qubits, for N even but not a multiple of 
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4. This is a system of N qubits on a circle with an interaction of the form XX + YY between 
nearest neighbour qubits, such that the overall Hamiltonian is 

'N-l 

Y, {XiX,+, + YiYi+i) + XiXn + YiYn 

. i=l 

The overall factor of j is included in order to match the original formulation in [18] and the 
notation used in the rest of literature. The properties of the system that we are most interested 
in are the spectral gap (the energy difference between the ground state |n) and the first excited 
state) and the spin correlation functions p = {Vt\XiXj + YiYj\Q) (the expectation value of 
XX + YY in the ground state). In particular we will show that both of these quantities are 
H(l/poly(iV)). 

This model has been extensively studied since the seminal paper by Lieb, Schultz and 
Mattis [18]. For a helpful review of the area see [21]. However, works in this area usually 
consider the infinite chain, whereas we will need to take N finite (but growing). In particular, 
asymptotic expressions for the spin correlation functions of the infinite chain were determined 
in [19]. Other related work studying correlation functions in varying settings (finite temperature, 
anisotropic interactions, etc.) includes [2, 20]. These can be approached using a mapping to 
the 2D Ising model (see e.g. [25]). Here we will start by considering the infinite chain, and 
then show that the finite chain provides an asymptotically precise approximation. Following 
the exposition in [21], we provide self-contained proofs of the results that we need. 



5.1 Ground State and Spectral Gap 

As is shown in [18] the Hamiltonian H can be transformed to a system of free fermions with 
nearest neighbour interactions in the following way: 

First let ^ ^ 

Qj = -{Xj + iYj) and = -{Xj - iYj) for j = 1,..., n 

so that the Hamiltonian of the system is now: 



'N-l 

^(ajoi+i -bOio|+i) 
_ i=\ 


+ ajojv T fflifljv 


Note that the ai do not satisfy canonical commutation relations for a boson or canonical 
anticommutation relations of a fermion, but instead a mix of the two: {ai, a\} = I, {ai, a,} = 0 
and [al,aj] = 0 = [ai,aj] for i j ■ So we apply a Jordan-Wigner transformation which takes 
the Hamiltonian into the form of one for a system of fermions, which can then be solved simply. 
Let 

/ *-i 

Ci = exp I iir 

V i=i 


and 


c| = a\ exp 


— ITT 



which satisfy the canonical fermionic anticommutation relations {ci,Cj} = 6ij and {ci,Cj{ = 
0. For more information on these operators see for example [21]. In terms of the c^’s the 
Hamiltonian becomes 



N-l N 

Y. + Cic\+i) - (cjvCi -h c\cn) exp(z7r ^ cjcj) 

i=l j=l 
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At this point Lieb, Schuitz and Mattis [18] turn to considering the simiiar Hamiitonian 

'N-l 

(cjci+l + Cic\+l) + cjvCl + c\^N 

. 

as they are interested oniy in the ieading order asymptotics, and this change oniy resuits in an 
0(1/A^) error in the caicuiation of the eigenvaiues. However, we wiii need to be more carefui 
here. 

Let P = exp(i7rAA) where J\f = fermionic number operator which has 

eigenvaiues 0,1, so the eigenvaiues of P are ±1. Note that P commutes with H, so we 

can diagonaiise H into biock diagonai form and consider its action on the two eigenspaces of P 
separateiy. Then 



H = H --— 


H< = 


'N-l 


^(c|ci+i + 


2=1 


where H- is as above, and 

Cicl+i) ~ (cjvCi + ({cn) 


Now we can diagonaiise both H- and separateiy, by first writing H- as ^ and 

finding the unitary matrix U that diagonaiises A. 




=► Aij — 


(h2j+5jAr)/2 i = \ 

{5i+ij + 6ji-i)/2 l<i<N 
{6ij + 6j i = N 

Let (jikj denote the jth coordinate of the kih. eigenvector of A. We can easiiy verify that with 
(pkj as beiow, {4>k\ is a compiete orthonormai set of eigenvectors with corresponding eigenvaiues 
Afc = cos(27r/c/Al) : 

1 


4’kj 


Vn 


^2'Kikj /N 


where 0 C k ^ N — I 


Then by making a second change of variabies rjk = Yli converseiy Cj = Ylk 4>ki'nk), 

which are aiso fermionic operators, the Hamiitonian becomes: 

H- = 

k 

To make the ground state energy ciear, we wiii make one finai change of variabies 

_ J r/fc Afc ^ 0 

^ I Afc < 0 

Then the Hamiitonian H- becomes 


^ ^k^Uk + ^ ^ lAfciCfcCfc +1 Afc I / 

k:Ak^0 k\Ak<0 k y/i::Afc<0 J 

And so the ground state of is |r2), the unique state that satisfies ^k\^) =0 for all fc, and 
which has energy equal to 


3N-2 

4 


fc:Afe<0 


COS 


n=^ 




r f 1 

/ 27rn\ 

= Re 

27ri(iV + 2) ^^ 27rin 

(a') 

e 4iv > e N 


n=0 
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= Re 


2 


1 


sin(f) 

The other eigenstates of H- will be of the form with extra energy above the 

ground energy of |Afc^| + + ... lA^J. Therefore in the case where N is not a multiple of 4, 

Afc / 0 for any k, and so |n) is the unique ground state of H-. 

In order for jf!) to also be a true eigenstate of the full Hamiltonian H with the same energy, 
we will also need P\0,) = —|H). Note that, for N even \{k : Ak < 0}| = N/2, so 

M\Q) = Y,cUQ) = Y,vivk\n)= \^) = j\^) 

i k k:Ak<0 

For N even and not a multiple of 4, N/2 is odd and so P\^) = exp(i7rAA)|n) = —|H) as required. 
We now repeat this procedure for the = '^clPijCj part of the Hamiltonian. 

( {S2j—^jN)/‘^ i = l 

H+ = '^ c\BijCj ^ Bij = < ((5j+i j + Sji-i)/2 1 <i < N 

( {-dij + 5jN-i)/2 i = N 

Let the eigenvalues and eigenvectors of B be given by A'^ and ^/kj respectively. They have 
a very similar form to A, only shifted by a factor of n/N: 

A/ = cos ('ipkj = , Oi^k^N-1 

V A' / pN 

In the case where N is even but not divisible by 4, the ground state of is fourfold degenerate 
and has energy 


leN 


27ri 

1 — e iv 




cos 


fc:A^sS0 


(2n + l)7r 
N 


= Re 




'-1 


E 

n=0 


2'!rin 

e N 


= Re 


I 27ri , 

L 1 — e v J 


tan(^) 


Since — sin(7r/A^) ^ < — tan(7r/A^) we have now shown that jH) is the unique ground 
state of the whole Hamiltonian H. 


5.1.1 Spectral gap 

In the P = — 1 sector, all states orthogonal to 
energy of at least | cos( ^*-^^^ )| = sin(^). 

And in the P = +1 sector, all states have 
least 

1 _ l__ ^ l-cos(^) 

sin(^) tan(^) sin (f) 

In particular one of the four eigenstates of P+ that has energy — tan(7r/A^)“^ will have 
P = +1, and so will be a true eigenstate of the total Hamiltonian H. Therefore the spectral 
gap is exactly tan (^). 


the |n) have an extra energy above the ground 
an extra energy above the ground energy of at 

2 sin (^) cos (^) V2Ary 
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5.2 Spin correlation functions 

We need to calculate the quantity {VL\XiXj + YiYj\Q), but given the symmetry of the problem 
it will suffice to calculate pij = {Q\XiXj\Q). Note that because of the translational invariance 
of the original Hamiltonian, we expect ptj to depend only on \i — j\. 

For arbitrary i < j, 


XiXj = (a| + ai){a] + aj) = (c| + a) exp^in ^ clck){c] + cj) 

k=i 


- ... - fJ 


= (cI-Ci) n ((4 +Cfc)(cl, -Cfc)) ( 

k=i-\-l 


c] + Cj 


So 


Pij — {^\BiAi+iBi+iAi+2 ■ ■ ■ Aj-iBj-iAj\Q) 


where we have defined the operators Ai = c| + = c| — Cj, which obey the following 

commutation rules: {Ai, Bj} = 0, {Ai, Aj} = 5ij, {Bi, Bj} = —25ij. We can use Wick’s Theorem 
to express this expectation in terms of a sum over all possible contractions. 

We have 


{Q\BiAm\fl) = (f^|(c| - Ci){cl, + C,„)|f^) = - (l)klVk){4>l'mVl + 

k,k' 


k y/c,A/5.<0 fc,A^^0 

Noting that this expression depends only on I — m, we define Gi-m = in order to 

match the formulation in [18]. Letting r = I — m and pairing up terms k and N — k, we get 

N-2 

4 2 

(i-(-ir)+2^ cos{2Trkr/N) — 2 cos{27rkr / N) 

k=l 

4 J 

Then we pair up the terms k and N/2 — k, which cancel out for even r, giving G,. = 0 if r is 
even. But for odd r: 



N-2 N ^ 

cos{2-Kkr/N) = - Y cos{2-Kkr/N) = -x + ^ 

k^ ^ 2sm(^) 


So 


/ V ■''+1 

_ 2(-l)^ 
iVsin(^) 


r odd 


Then from Wick’s Theorem (and the calculations {AiAm) = ^im = {BiBm) which can be 
found in [18]), we get pij = det Rn,N, where n = i — j and Rn,N is the nxn matrix with entries: 


{Rn,N)lm = Gl —m+1 


(3) 
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In [18] (and the rest of the literature) the limit —>■ oo is taken earlier as this simplifies the 
calculation of Gi-m (especially in the anisotropic case), and again does not affect the leading 
order asymptotics. This then gives a matrix Rn, with entries: 


(Rn) 


Im — 


. . l — m 

(- 1 ) — 


7r(/—m+1) 


l — m odd 
I — m even 


(4) 


Given the known limit lima;_^o x / sin x = 1 we see that each entry in Rn is the limit as 
—)• oo of the corresponding entry in Rn,N- Since det is a polynomial, and hence continuous, 
function of the entries of a matrix, for fixed n 

det Rn = lim deti?„Ar. 

N^oo ’ 

However we will need a stronger result for the behaviour as n —)• oo and N = poly(n). The exact 
result is given in Lemma 8 below. We will first need to determine the determinant of the matrix 
Rn- This was previously obtained in [19]; for completeness, we include an alternative concise 
proof based on the general theory of Toeplitz matrices and a recent result of Ehrhardt [14] on 
the Fisher-Hartwig conjecture. 


5.3 Toeplitz matrices 

. as a Toeplitz matrix, see for example [4] for 

^2 • ■ • tn ^ 

h : 

to 

to ti 

t-i to 

These matrices have been especially well studied when there exists a complex-valued function 
/ integrable on the unit circle such that {t^} are the Fourier coefficients of /: 

4 = ^ / kGZ (5) 

Jo 

The function / is called the symbol of the corresponding Toeplitz matrix Rn- We can find such 
a symbol by calculating the Fourier series with coefficients t^- If f{z) is sufficiently smooth, 
then the famous Szego limit theorem gives the asymptotic behaviour of det(i?n) as n —>■ oo. 
However in our case, the symbol of Rn is given by the function g{z), which is discontinuous: 

oo ( z arg(z) G [0, |) U (^, 27r) 

aiz) = tkz’^ = < -z arg( 2 :) G (f, ^) (6) 

k=—oo V ® 2: = zbi 

Note that this function g is only defined on the unit circle S = {z \z\ = 1}. In fact the value 
of g{z) at z = is unimportant as the coefficients tk as defined by (5) will be the same for 
functions that agree almost everywhere (i.e. functions that differ only on a set of Lebesgue 
measure zero). 

The Fisher-Hartwig conjecture is an attempt to generalise the Szego limit theorem to piece- 
wise continuous symbols, such as g{z)] see [4] and references therein. The Fisher-Hartwig 


A matrix of the following diagonal form is know] 
an introduction. 

/ to ti 


Tn = 


t-l to 
t-2 t-l 

\ t-n --- 
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conjecture has previously been used to study correlation functions of the XY model, for exam¬ 
ple in [23]. A function / is said to have Fisher-Hartwig singularities if it can be written in the 
following form: 


f{z) = Yl\z- z,\^'-^g,^,p^{z)z/^ 

j=0 


(7) 


Q ^ arg(2;) < 9j 
g-* 7 r/ 3 j 0 ^ ^ arg(2;) < 27r 


where zj = e*®-’ , 9j G [0,27r), and gzj,i3j = 

The Fisher-Hartwig conjecture states that for V{z) sufficiently smooth and otj,^j G C, then 
det(i?„) ~ E[f] exp{nVo)n^i^o as n —)■ oo 

where E[f] is a constant independent of n, and Vo = ^ j‘^^V{E^)d9 is the zero-th Fourier 
coefficient of V. 

This conjecture has now been proved to be true for a number of different conditions on aj, Pj 
and V; but in particular Ehrhardt [14] showed that this holds under the following conditions: 

i) Re aj > —1/2 for all j 

ii) maxj^k |Re /3j — Re /3fcj < 1 for all j, k 

iii) aj ± f3j / —1, —2, —3,... for all j 

iv) V G C°° {V is infinitely differentiable). 


An alternative proof of this result is presented in [13] . 

The function g{z) can be written in the form of (7) with two Fisher-Hartwig singularities: 

g{z) = zgi_i{z)g_-_i{z) 

Note that for j = 0, 1, aj = 0,l3j = —1/2 and H = 0, so all of the four conditions (i)-(iv) are 
satisfied. In fact the exact constant can also be calculated: 

det(i?ri) ~ En~2 

where E = ^/2G{\fG{lf = ^ « 0.58835... 

where G is the Barnes G-function and A is the Glaisher-Kinkelin constant. This agrees with 
the result previously obtained in [19]. 


5.4 Relation between Rn,N and 

We now prove the following Lemma: 

Lemma 8. Let Rn,N o-nd Rn he the n x n matrices defined in equations (3) and (4)- Then, if 
N = Vt{n^) for a > 7/4, 

det(R„^ 7 v) ~ det(i2n) as n ^ oo. 

Remark. We have strong numerical evidence to suggest that this result still holds under 
more relaxed conditions for N, perhaps even just N > n. However the result as stated here is 
enough for our purposes. 
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Proof. The proof uses the standard Weyl’s inequality for Hermitian matrices to relate the 
eigenvalues of matrices closely related to Rn and Rn,N- First introduce the permutation matrix 
Sn with entries {Sn)ij = 5i^n+i-j- Then RnSn and Rn^NSn are both real symmetric matrices. 
Let {Aj} and {fii} be the ordered eigenvalues of RnSn and Rn^NSn respectively. Then, by Weyl’s 
inequality: 

~ ll-®ll ^ ^ + 11-611 

where B = Rn,N — Rn is the difference between these two matrices. Then divide through by A* 
(considering the cases A* > 0 and Aj < 0 separately), to obtain 


1 - 


\\B\ 




1 + 


\\B\ 


From the general theory of Toeplitz matrices with bounded symbols, we know that ||i?n|| ^ 
Halloo = 1, where g is as defined in (6). So for all i, |Aj| ^ ll-^nS'll = ||.Rn|| ^ 1- Also, 
Hi Aj = det(i2ri) as re —>■ oo, so there exists E' such that for sufficiently large re, 

|Aj| > E'. Then taking the product over v. 


1 - E'-^\\B\\n 


1/2 


" / A Aij det(i?n,Ar) 


n 

2=1 




Xi det(2?,2) 


^ il + E'-^\\B\\n 


A/2 


( 8 ) 


Note that B is another Toeplitz matrix, with entries Bij = where 6^ = 0 for even r, 

and for odd r : 


hr = (-1) 2 — 


o / Trr 

1 2 / ^ 


- 1 


r odd 


vrr Vsin(^) 

To deal with the term in brackets, we define h{x) = x/ sinx — 1. Note that since /i^(0) = 
0, /i"(0) = 1/3 and h"{x) is continuous in a region around zero, there exists C > 0 and (i > 0 
such that for all x G {—5,5)'. 

2 

h{x) = —— 1 = — + 0{x^) < Cx^. 
sm(x) 6 

Then for all i and j, and N ^ TTnj5, we can bound \Bij\ ^ |6„| < 2CTm/N‘^ so we can upper 
bound the operator norm ||i?||oo, by the Hilbert-Schmidt norm: 

2Ctt 


^ WBWo = 


E i%i" < 


\ ij-i 


iV2 


, E 

\ *,i=i 


re^ = 


2C7rre^ 
iV2 ' 


So if = fl(re“), for some a > 7/4, then there exists a constant D > 0 such that 


/ \ '^ / JyTl, 2 \ ^ 

^1 + ^ [ 1 H-I ^ exp(T>re 2 “^“) —> 1 as re —> oo 

where we have used the inequality 1 + x ^ e^. Similarly, we can use the inequality 1 — x ^ 
which is valid for x G [0, e] for some e > 0. Then for re sufficiently large 


1 - E'-^WBWn^^A"' A ( 1 - 


Dn2 


-2a 


n 


A exp(— 2 iAre 2 1 as re —)> oo 


Finally, combining these results into equation (8), we get 

det(i?„,Ar) 
lim ——^—E— = 1. 
n-^oo det(it„) 

□ 
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The following lemma summarises the above discussion: 

Lemma 6 (restated). Fix N even but not a multiple of 4, o.nd let H = + 

YiYij^i)+XiXN+YiYi^. Then H has a nondegenerate ground state |n) and spectral gap Tl{l/N). 
Further, for any pair i, j such that \i—j\ = n andn = o{N‘^/'^), {Q\XiXj + YiYj\il.) = 

There is an efficient classical algorithm to compute the spectral gap and all the correlation 
functions. 

6 Outlook 

Although we have translated some previous results about QMA-completeness of the Local 
Hamiltonian problem closer to truly physically realistic systems, one significant issue remain¬ 
ing is the large weights required for the interactions in the hard instances. This requirement 
is a basic limitation imposed by our use of perturbation theory for gadgets. However, it would 
be very interesting if recently developed techniques proving QMA-hardness with lower-weight 
interactions [9, 10, 11] could be extended or combined with our results. 

There are also some specific open questions left over. First, the natural cases of the general 
Heisenberg model on a square lattice, and the antiferromagnetic Heisenberg model on a trian¬ 
gular lattice, are still unresolved. Proving these cases QMA-complete could follow from finding 
an exactly solvable (but nontrivial) special case of the Heisenberg model whose interactions are 
suitably sparse. Second, there is scope for tightening our classification of interactions of the 
form aXX -|- fiYY -|- jZZ with fixed signs. For many of the interactions of this form, we know 
the interaction is contained within StoqMA, but not whether it is StoqMA-complete, or within 
some smaller complexity class. A particularly natural example which is currently unknown is 
the ferromagnetic XY model, whose interactions are of the form —XX — YY. Third, we have 
not been able to resolve the complexity of every set of 2-qubit interactions with fixed signs. It 
remains to classify interactions with nonzero 1-local part, and the case of the 5'’'-Hamiltonian 
problem where |5| ^ 2. 

These open questions, along with our results described here, highlight that the Local 
Hamiltonian problem displays a rich and complex structure when restrictions on interaction 
signs and topologies are considered. 
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